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Abstract 

Two computational models to determine the fatigue life and reliability of a commercial turboprop 
gearbox are compared with each other and with field data. These models are (1) Monte Carlo simulation 
of randomly selected lives of individual bearings and gears comprising the system and (2) two-parameter 
Weibull distribution function for bearings and gears comprising the system using strict-series system 
reliability to combine the calculated individual component lives in the gearbox. The Monte Carlo 
simulation included the virtual testing of 744,450 gearboxes. Two sets of field data were obtained from 64 
gearboxes that were first-run to removal for cause, were refurbished and placed back in service, and then 
were second-run until removal for cause. A series of equations were empirically developed from the 
Monte Carlo simulation to determine the statistical variation in predicted life and Weibull slope as a 
function of the number of gearboxes failed. The resultant L w life from the field data was 5,627 hr. From 
strict-series system reliability, the predicted L 10 life was 11 A hr. From the Monte Carlo simulation, the 
median value for the L 10 gearbox lives equaled 757 hr. Flalf of the gearbox L 10 lives will be less than this 
value and the other half more. The resultant L UI life of the second-run (refurbished) gearboxes was 
1,334 hr. The apparent load-life exponent p for the roller bearings is 5.2. Were the bearing lives to be 
recalculated with a load-life exponent p equal to 5.2, the predicted L w life of the gearbox would be equal 
to the actual life obtained in the field. The component failure distribution of the gearbox from the Monte 
Carlo simulation was nearly identical to that using the strict-series system reliability analysis, proving the 
compatibility of these methods. 


Nomenclature 


ci i life factor for reliability 
a 2 life factor for material and processing 

a 2 life factor for operating conditions (lubrication and environment) 
B material life factor or constant for gear life 
C D dynamic load capacity, N (lbf) 
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C, basic load capacity of gear tooth, N (lbf) 
c stress-life exponent, N (lbf) 

d diameter of rolling element, mm (in.) 
e Weibull slope 

/ face width of gear tooth, m (in.) 

f cm bearing geometry and material coefficient 

h exponent 

k number of load cycles of gear tooth per input shaft revolution 

L life, number of stress cycles or hr 
L„ component L 10 life 

Lp characteristic life or life at which 63.2 percent of population fails, number of stress cycles or hr 

L 10 10-percent life or life at which 90 percent of a population survives, number of stress cycles or hr 

L 5 o 50-percent life or life at which 50 percent of a population survives, number of stress cycles or hr 

/ involute length, m (in.) 

m Weibull modulus (same as Weibull slope e) 

N number of gear teeth on a gear 

n number of components or gearboxes failed 
P eq equivalent radial load, N (lbf) 

P, normal tooth load, N (lbf) 

p load-life exponent 

r\ pitch radius of pinion, m (in.) 

r 2 pitch radius of gear, m (in.) 

S probability of survival, fraction or percent 
V stressed volume, m 3 (in. 3 ) 

X fractional percent of components failed 

X„ fractional percent of time at operational condition 

Z number of rolling elements 

z depth below surface, m (in.) 

z 0 depth to orthogonal shearing stress, m (in.) 

a contact angle, deg 

r| number of stress cycles to failure 

r| 10 life for single gear tooth at 90-percent probability of survival, millions of stress cycles 

p curvature sum, m-1 (in.-l) = (1/ri + l/r 2 )(l/sincp) 
x 0 orthogonal shearing stress, GPa (ksi) 

cp pressure angle, deg 

Subscripts: 

B bearing 

G gear 

i i* component 

sys system probability of survival or life 

t gear tooth 

P designation of characteristic life 

1,2 bodies 1 or 2; load-life condition 1, 2, etc. 
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Introduction 


Predicting and verifying the life of aircraft mechanical power trains, such as turboprop and helicopter 
main rotor gearboxes, is a complex task. Design requirements such as low weight and high power 
capacity need to be balanced with long life, high reliability, and low maintenance costs. The primary 
components of these power trains are combinations of rolling-element bearings, gears, and shafting. 
Accurate prediction of the lives of these components is necessary to establish reliability, predict 
replacement rates, set up maintenance and replacement schedules, and ascertain warranty limits. 

The lives of these mechanical components are not deterministic but are probabilistic. A significant 
cause for power train removal and repair is contact (rolling-element) surface fatigue of the bearings 
and/or of the gears. Rolling-element fatigue is extremely variable but is statistically predictable depending 
on material (steel) type, processing, manufacturing, and operating conditions. 

Rolling-element fatigue life analysis is based on the initiation or first evidence of fatigue spalling on a 
loaded, contacting surface of either a bearing or a gear. This spalling phenomenon is load cycle 
dependent. Generally, the spall begins in the region of maximum shear stresses, located below the contact 
surface, and propagates into a crack network. Failures other than that caused by classical rolling-element 
fatigue are considered avoidable if the component is properly designed, handled, installed, and not 
overloaded (ref. 1). Flowever, under low elastohydrodynamic (EHD) lubricant film conditions, rolling- 
element fatigue can be surface or near-surface initiated with the spall propagating into the region of 
maximum shearing stresses. 

In 1939, Weibull (refs. 2 to 4) developed a method and an equation for statistically evaluating the 
fracture strength of materials based upon small sample (population) sizes. This method can be and has 
been applied to analyze, determine, and predict the cumulative statistical distribution of fatigue failure or 
any other phenomenon or physical characteristic that manifests a statistical distribution. 

Based on the work of Weibull (ref. 2) in 1939, Lundberg and Palmgren (ref. 5) in 1947 showed that 
the probability of survival S could be expressed as a power function of the orthogonal shear stress r c , life 
r\ e , depth to the maximum orthogonal shear stress z 0 , and stressed volume V. That is, 


ln- 


-V 


From equation (1), Lundberg and Palmgren (ref. 5) formulated the following relation for rolling- 
element bearings: 


MO 





where Co is the dynamic load capacity (the load that a bearing can carry for a life of 1 million inner-race 
revolutions with a 90-percent probability of survival (L 10 life)), P eq is the equivalent radial load, and p is 
the load life exponent. 

Where the data base for gear life is limited, the data base for ball and roller bearings is significantly 
more extensive. A concern that arises from these data and the analysis is the variation between life 
calculations and the actual endurance characteristics of these components. Experience has shown that 
endurance tests of groups of identical bearings under identical conditions can produce a variation in L\ 0 
life from group to group. If a number of apparently identical bearings are tested to fatigue at a specific 
load, there is a wide dispersion of life among these bearings. For a group of 30 or more bearings, the ratio 
of the longest to the shortest life may be 20 or more (ref. 1). For gears there is a similar variability. This 
variation can exceed reasonable engineering expectations. 
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Between 1975 and 1981, Coy, Townsend, and Zaretsky (refs. 6 to 9) published a series of papers 
developing a methodology for calculating the life of spur and helical gears based upon the Lundberg- 
Palmgren theory and methodology for rolling-element bearings. These gear-life models utilize the 
Lundberg-Palmgren theory. Also developed from these models is a reliability model for planetary gear 
trains (refs. 10 to 15), which takes into account the individual reliabilities of the gearbox bearing and 
gears based upon classical rolling-element fatigue. The reliability of the gearbox system is treated as a 
strict-series probability combination of the reliabilities of the gearbox components in accordance with 
Lundberg-Palmgren theory (ref. 5). Each bearing and gear life was calculated and the results statistically 
combined to produce a system life for the total gearbox. Because of the different component Weibull 
slopes, a linear approximation was made to determine the system Weibull slope. The method was applied 
to a turboprop gearbox by Lewicki et al. (ref. 16). 

Another method for determining the life of a system comprising multiple components is the Monte 
Carlo simulation (ref. 17). Vlcek et al. (ref. 18) determined the life of rolling-element bearings using this 
method and the two-parameter Weibull distribution function (ref. 4). 

In view of the aforementioned, the following are objectives of the work reported herein: (a) to 
determine the variation in gearbox lives and distribution parameters as a function of sample size using a 
Monte Carlo simulation of randomly selected lives of individual bearings and gears comprising the 
system, (b) to determine the life of a gearbox based on a two-parameter Weibull distribution function for 
bearings and gears comprising the system using strict-series system reliability, and (c) to compare both 
methods with each other and with component and gearbox lives obtained from field data. 


Procedure and Analysis 

Turboprop Gearbox 

The commercial turboprop gearbox used for this analysis (fig. 1) consists of two stages with a single- 
mesh spur reduction followed by a 5-planet planetary gearbox consisting of 1 1 rolling-element bearings 
and 9 spur gears (ref. 16). The first stage consists of the input pinion gear meshing with the main drive 
gear. The second stage is provided by the fixed-ring planetary driven by a floating sun gear as input with 
a five-planet carrier as output. At cruise conditions, the input pinion speed is constant at 13,820 lpm, 
producing a carrier output speed of 1 ,02 1 rpm. A list of the component parts of the gearbox is given in 
table 1. 


Front pinion 
roller bearing 

Ring gear-. 
Planet spherical 
roller bearing — 

Planet carrier—,. 

Planet gear—. 


Output shaft 

Radial cylindrical- 
bearing 

Thrust ball bearing- 



Pinion gear 
,^Rear pinion 
roller bearing 

Input shaft 

Main drive gear 

— Main drive 
roller bearing 

Carrier support 
roller bearing 


'-Sun gearl 


Figure 1. — Commercial turboprop gearbox. 


NASA/TM— 2007-215019 


4 



TABLE 1.— PREDICTED TURBOPROP GEARBOX COMPONENT LIVES FROM 
LUNDBERG-PALMGREN ANALYSIS WITH LIFE FACTORS AND 
STRICT-SERIES SYSTEM RELIABILITY 


Component description 

| Predicted life, hr 

Weibull 

(see fig. 1 ) 

1 L\ o 1 

C 50 

slope, e 

| Rolling-element bearings | 

Cylindrical roller bearings 
Front pinion 

20,890 

111,476 

1 . 

125 

Rear pinion 

21,312 

113,728 



Main drive 

26,459 

141,194 



Carrier support 

312,881 

1,669,635 



Propeller radial 

68,194 

363,905 


r 

Propeller thrust ball bearing 

33,065 

180,484 

1. 

11 

Planet double-row spherical roller bearing set, 5 bearings 

844 

4,503 

1.125 

Bearing system life 

774 

4,132 

1.125 

| Gears ! 

Pinion 

53,477 

131,552 

1 2 - 5 I 

Ring 

4,540,212 

11,168,760 



Sun 

19,033 

46,821 



Main drive 

108,148 

266,040 



Planet gear set, 5 gears 

28,092 

69,105 

1 

r 

Gear system life 

16,680 

44,032 

1 2.5 i 

Gearbox life 

774 

4,132 

1 1. 

.125 


A typical mission profile for this commercial gearbox is given in table 2, which presents the duration 
as a percentage and the propeller shaft power for each flight condition. This profile includes loads for (a) 
takeoff, (b) climb, (c) cruise, and (d) descent. The cruise segment of the profile consumes 68 percent of 
the flight time with a little less than half of the power required for the takeoff, which lasts for less than 
3 percent of the flight time. 


TABLE 2.— MISSION PROFILE OF COMMERCIAL 


TURBOPROP GEARBOX 


Mission 

segment 

Percent time 
of segment 

Propeller shaft power, 
kW (hp) 

Takeoff 

2.84 

3,132 (4,200) 

Climb 

17.02 

2,461 (3,300) 

Cruise 

68.08 

1,516(2,033) 

Descent 

12.06 

945 (1,267) 

Equivalent 

100.00 

1,833 (2,457) 


The cause for removal can be assumed to be one or more bearings or gears that had fatigue or damage 
resulting in wear and/or vibration detected by magnetic chip detectors and/or vibration pickups. The 
gearbox is removed from service before secondary damage occurs and is inspected. After the failed part 
or parts are replaced, the gearbox is put back into service. 

Individual failure occurrences are not predictable but are probabilistic. No two gearboxes run under 
the same conditions fail necessarily from the same cause and/or at the same time. At a given probability 
of survival, the life of the gearbox system will always be less than that of the lowest lived element in it. 

Historical field data for 64 gearboxes were collected. The first part of these data covered the time 
from their field installation and first field operation to their removal for cause (failure) and refurbishment. 
The second part covered the time after refurbishment and placement back into field operation of 59 of the 
64 gearboxes until their second removal for cause. 
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Figure 2. — Weibull plot where (Weibull) slope 
of tangent of line is e; probability of survival, 

Sp, is 36.8 percent at which L = Lp, or L/Lp = 1 . 


Weibull Analysis 

In 1939, Weibull (refs. 2 to 4) developed a method and an equation for statistically evaluating the 
fracture strength of materials. He also applied the method and equation to fatigue data based upon small 
sample (population) sizes, where the two-parameter expression relating life and probability of survival is 


f i V 


lnln 


= eln 


f l' 






where 0<L<co;0<5'<1 


( 3 ) 


When plotting the In ln( IAS’) as the ordinate against the In L as the abscissa, fatigue data are assumed to 
plot as a straight line (fig. 2). The ordinate In ln( IAS’) is graduated in statistical percent of components 
failed or removed for cause as a function of In L, the natural logarithm of the time or cycles to failure. The 
tangent of the line is designated thOe Weibull slope e, which is indicative of the shape of the cumulative 
distribution or the amount of scatter of the data. The method of using the Weibull distribution function for 
data analysis to determine component life and reliability was later developed and refined by L.G. Johnson 
(ref. 19). This method was used to analyze the data reported herein. 


Bearing Life Analysis 

Lundberg and Palmgren (refs. 5 and 20) extended the theoretical work of Weibull (refs. 2 to 4) and 
showed that the probability of survival S could be expressed as a power function of shear stress x, number 
of stress cycles q, depth of maximum shear stress z, and stressed volume V (see eq. (1)). By substituting 
the bearing geometry and the Hertzian contact stresses for a given load into equation (3), the bearing basic 
dynamic load capacity C D can be calculated (ref. 5). Lundberg and Palmgren related the L I0 life to the 
equivalent bearing load, load-life exponent, and the basic dynamic load capacity Co (eq. (2)). 

For ball bearings, the generally accepted value for the exponent p in equation (2) is 3 (ref. 5). 

Initially, for roller bearings the derived value for p was 4 (ref. 5). However, in their subsequent work, 
Lundberg and Palmgren (ref. 20) recommended a value of p of 10/3. The ANSI/ABMA and ISO 
standards (refs. 21 to 23) adopted a value of p equal to 10/3 for roller bearings and a value of 3 for ball 
bearings. Bearing computer codes in use today use p = 3 for ball bearings and p = 4 for roller bearings. 
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Formulas for the basic dynamic load ratings derived by Lundberg and Palmgren (refs. 5 and 20) are 
incorporated in the ANSI/AFBMA and ISO standards (refs. 21 and 23). 

Equation (2) can be modified using life factors based on reliability a\, materials and processing a 2 , 
and operating conditions such as lubrication a 2 (ref. 1) where 

h — 


Gear Life Analysis 

Between 1975 and 1981, Coy, Townsend, and Zaretsky (refs. 7 to 9) published a series of papers 
developing a methodology for calculating the life of spur and helical gears based upon the Lundberg- 
Palmgren theory and methodology for rolling-element bearings. Townsend, Coy, and Zaretsky (ref. 24) 
reported that for AISI 9310 spur gears, the Weibull slope e is 2.5. Based on equation (1), for all gears 
except a planet gear, the gear life can be written as 

* 


For a planet gear, the life is 


Ao- — ' 


N 


~ l/e G i-n~ e G 


Tio,, +P 


~ e G 

10 ,-, 


|-!/e G 


r 


k 


The A | o life of a single gear tooth can be written as 


bio, 


: a 2 a 3 


r c^ PG 

vAy 


where 


/ i D r0.901 - 1 . 165 ,- 0.093 

C t =Bf p / 


and 


P = 


1 1 



1 

sin cp 


and r)io, is the L UI life in millions of stress cycles for one particular gear tooth. This number can be 

determined by using equation (7) where for bearings, C, is the basic load capacity of the gear tooth; P, is 
the normal tooth load; p G is the load-life exponent usually taken as 4.3 for gears based on experimental 
data for AISI 9310 steel; and d 2 and a 2 are life adjustment factors similar to those for rolling-element 
bearings. The value for C t can be determined by using equation (8), where B is a material constant that is 
based on experimental data and is approximately equal to 1.39xl0 8 when calculating C, in SI units 
(Newtons and meters) and is 21,800 in English units (pounds and inches) for AISI 9310 steel spur gears;/ 
is the tooth width; and p is the curvature sum at the start of single-tooth contact. 

The / 0f life of the gear (all teeth) in millions of input shaft revolutions at which 90 percent will survive 

can be determined from equation (5) or (6) where A is the total number of teeth on the gear, e ( , is the 
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Weibull slope for the gear and is assumed to be 2.5 (from (ref. 24)), and k is the number of load (stress) 
cycles on a gear tooth per input shaft revolution. 

For all gears except the planet gears, each tooth will see a load on only one side of its face for a given 
direction of input shaft rotation. Flowever, each tooth on a planet gear will see contact on both sides of its 
face for a given direction of input shaft rotation. One side of its face will contact a tooth on the sun gear, 
and the other side of its face will contact a tooth on the ring gear. Equation (6) takes this into account, 
where r| 10 is the L ]() life in millions of stress cycles of a planet tooth meshing with the sun gear, and 

r| 10(7 is the T 10 life in millions of stress cycles of a planet tooth meshing with the ring gear. 

System Life Prediction 

The L io lives of the individual bearings and gears that make up a rotating machine are calculated for 
each condition of their operating profiles. For each component, the resulting lives from each of the 
operating conditions are combined using the linear damage (Palmgren-Langer-Miner) rule (refs. 25 to 27) 
where 


L L x L 2 L n 


( 10 ) 


The lives of each of the gearbox components are combined to determine the calculated system L l0 life 
using the two-parameter Weibull distribution function (eq. (3)) for the bearings and gears comprising the 
system and strict-series system reliability (ref. 10) as follows: 


— -Y — - 

*e Z-f re 


1 1 
+ + • 


sys 


i = 1 


r t 


n J 


(lla) 


The system shape parameter (Weibull slope) using strict-series reliability (eq. (lla)) where each 
component or combinations of multiple bearings and gears have different Weibull slopes (shape 
parameters) is not intuitively obvious. If at each time sequence the probability of survival of each 
component is multiplied together, the system reliability at that time and, hence, the probability of failure 
can be determined. When these values are plotted on Weibull paper, a Weibull slope (shape parameter) 
can be determined for the system life distribution using a least squares fit. When this is done, it is found 
that the system Weibull slope approximates that of the lowest lived component in the system. For the 
gearbox, 


f 


sys 


T e B T e B 

D D 


+ ...- 


V 


L e * 

a n J 


\ r 

1 + 


1 1 

T e G T e G 

V C, G 2 


+ ...- 


L ? 

G„ j 


(Hb) 


The lowest lived components in the gearbox are the roller bearings. As a result, the Weibull slope 
assumed for the planetary gear spherical roller bearings is assumed to be the Weibull slope e of the entire 
gearbox system in eq. (lib). 
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Virtual Gearbox Testing 

It was assumed that the virtual gearbox was assembled from 12 virtual bins containing components 
listed in table 1. As in a realistic manufacturing process, each of the component parts of the gearbox is 
grouped in separate bins. Each of the virtual bins contains 1,000 of the respective components. According 
to the method of Johnson (ref. 19), each component part in the bin is assigned an order number (1,2, 

3, ... 1 ,000) which is also the order in which they would fail. The life of each component at the probability 
of survival corresponding to the randomly selected ra nk can be determined using equation (3), knowing 
the survivability (from the order and ranking), characteristic life, and Weibull slope. 

Using Monte Carlo techniques (ref. 18), various-sized groups or populations of gearboxes were 
randomly assembled from the virtual part bins. This group or population represented failure distributions 
of varying sizes; sizes of 2, 3, 4, 5, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 35, 40, 45, 50, 58, 72, 
100, 200, and 1000 were arbitrarily selected. At population sizes greater than 20, the interval between 
population size was increased. Each random number generated within a population corresponded to an 
order number of a component in the virtual bin, with which an associated failure life was determined. 

The life of the gearbox was that of the shortest lived component. All gearboxes in each set failed; 
there were no suspensions or censoring. As a result of the population sizes evaluated, a total of 744,450 
virtual gearboxes was assembled and failed in this manner. From each population size n, a distribution of 
associated lives was thus determined. 

Next, this distribution of fatigue lives was sequentially ranked and then ordered according to the 
method of Johnson (ref. 19). From the order number, a percent survivability was determined 
corresponding to each fatigue life. A curve fit of a Weibull plot with the independent variable of the 
natural logarithm of the life (In L) and a dependent variable of the natural logarithm of the natural 
logarithm of the inverse of the probability of survival (In ln( IAS')) was created. A straight line was fitted to 
these n data points using the method of least squares. The Weibull slope and the L\ t] life of the gearbox 
systems were determined from the resultant line. For each value of n, the procedure was repeated 21 times 
to estimate variation between trials and to determine the maximum, minimum, and median values for 
each series of 21 trials. The 1 1th trial in ascending order of the 21 trials represents the median value. 
Ninety percent of the trials (19 trials) were bounded by the maximum and minimum L l0 life values. This 
entire process in turn was repeated 10 times to establish statistical variation between repetitions. Thus, the 
lives and slopes reported herein are averaged values from these 10 repetitions. 

Results and Discussion 

Predicted Gearbox Life 

The commercial turboprop gearbox used for this analysis is shown in figure 1. Using equations (2) 
and (3) for bearings and equations (5) through (9) for gears and appropriate computer programs 
incorporating these equations, the lives of each of the bearings and gears making up the gearbox were 
predicted for each of their operating conditions. These lives are shown as the Weibull plots in figures 3(a) 
and (b), respectively. 

The L io life of a single double-row spherical bearing is 3,529 hr. From equation (11a), the system L\q 
life for the five-bearing planetary set is 11 A hr. For all the bearings in the gearbox, the bearing system L l0 
life is also predicted to be 774 hr. 

Using equation (17a) for the individual gears, the gear system predicted L w life is 16,680 hr. 
Combining the bearing and gear lives to obtain a gearbox L 10 life, using equation (17b), the predicated L Ul 
life for the gearbox is 11 A hr. The lives of the individual bearings, and more specifically, those of the 
planet double -row spherical bearings, determine the life of the gearbox in this example. The system lives 
of the bearings, gears, and gearbox are summarized in figure 3(c). 
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Statistical percent of component Statistical percent of gears failed Statistical percent of bearings failed 

system failed 


95 
90 
60 
40 

20 h 

10 
6 
4 


Planet double 
row spherical 
roller bearing 
set, 5 each 



Prop thrust 
ball bearing 


Bearing 
system 
life — 


(a) 


100 


1 
1,000 



Main drive 
cylindrical roller 
bearing 
Rear pinion 
cylindrical 
bearings 

Front pinion 
cylindrical roller 
bearing — - 


Carrier support 
cylindrical 
roller bearing 


10,000 100,000 

Time, hr 


Prop radial 
cylindrical 
roller bearing 

I 

1 , 000,000 


10 , 000,000 




Time, hr 


Figure 3. — Predicted lives for commercial turboprop gearbox and its respective bearing and gear 
components using Lundberg-Palmgren life theory, (a) Bearing component lives, (b) Gear component 
lives, (c) Gearbox life and component lives. 
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Virtual Gearbox Life 


For this study, gearboxes were randomly assembled by Monte Carlo (random) number generation 
(ref. 29) from 12 virtual bins. The method of Johnson (ref. 19) was used to determine the Weibull slope 
and lives. For each population size, the process to determine the virtual gearbox lives was repeated 21 
times, from which the maximum and minimum lives and slopes could be determined at each population 
size. In turn, the process was repeated 10 times to establish the statistical variation, and the virtual lives 
and slopes are reported as averages of the 10 repetitions. The maximum and minimum lives of 21 trials 
were averaged for 10 repetitions and plotted as a function of population size in figure 4(a). The median 
values for the L\ 0 gearbox lives equal 757 hr and reasonably correlate with the predicted value of 11 A hr 
from the Lundberg-Palmgren method. Half the gearbox Z 10 lives, however, will be less than this value and 
the other half greater. 




Figure 4. — Average values of maximum and minimum 
variations of life and Weibull slope from Monte Carlo 
simulation of assembly and testing of 744,450 com- 
mercial turboprop gearboxes as a function of number 
of gearboxes failed, (a) Go life, (b) Weibull slope. 
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The variation in life is a function of the number of gearboxes failed and not the number of gearboxes 
tested or in the field (ref. 28). The equations providing a reasonable bounding fit of these data (fig. 4) are 
as follows: 


Maximum L \ q life = (median or predicted L\q bfe)(l + 6.5 n 0 5 ) 

Minimum Z 40 life = (median or predicted L\q life)(l.O-1.85« -0 - 3 ) 

Minimum L\ 0 life — > 0 where n< 3 

The maximum and minimum Weibull slopes of 21 trials were also averaged for 10 repetitions and 
plotted as a function of population size in figure 4(b). The variations in Weibull slope as a function of the 
number of gearboxes failed are as follows: 

Maximum Weibull slope = 1 . 125 + [5.0(ln«)]~°' 5 
Minimum Weibull slope = 1 . 125 - ( l . 125 m ) -0 ' 3 


(12a) 

(12b) 

(12c) 


(13a) 

(13b) 


There is good correlation between the median values of the Weibull slopes from the Monte Carlo 
simulation and that predicted for the gearboxes based on the Weibull slope of the lowest lived component 
in the gearbox itself, particularly for trial populations greater than 58 (table 3). However, half the resultant 
Weibull slopes will be greater than that predicted and the other half less. Again, the magnitude of the 
Weibull slope variation depends on the number of gearboxes that fail in a set and not the number in a set 
(refs. 19 and 28). 


TABLE 3.— SUMMARY OF LIVES FROM FIELD DATA 
OF NEW AND REFURBISHED COMMERCIAL 
TURBOPROP GEARBOXES COMPARED WITH 
PREDICTED LIVES USING LUNBERG- 
PALMGREN LIFE THEORY 


Life, hr 

Field failure data gearbox 
designation 

Predicted 
(table 1) 

New 

Refurbished 

£10 

5,627 

1,334 

774 

L50 

13,303 

5,597 

4,132 

Characteristic, Lp 

15,727 

7,398 

5,721 

Mean 

a 13,928 

b 6,819 

c 5,48 1 

Weibull slope, e 

2.189 

1.313 

1.125 


d At a 53.5-percent probability of failure. 
b At a 59.3-percent probability of failure. 
c At a 61.4-percent probability of failure. 


The work reported herein shows that gearbox or drive train failure is probabilistic and not 
deterministic. The specific life of individual drive trains cannot be determined with any reasonable 
engineering certainty. Further, the remaining life of a drive train in service cannot be determined prior to 
actual failure of a component within the gearbox. Assuming an ideal and accurate gearbox life prediction 
procedure, half of all gearbox groups from carefully controlled assemblies will have lives less than that 
predicted. However, the variation in the resultant life from that predicted will decrease as the number of 
failed gearboxes in the randomly selected group increases. 
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Gearbox Field Data 


New gearboxes 

The application of the Lundberg-Palmgren theory (ref. 5) to predict gearbox life and reliability needs 
to be benchmarked and verified under a varied load and operating profile. The cost and time to laboratory 
test a statistically significant number of gearboxes to determine their life and reliability is prohibitive. A 
practical solution to this problem is to benchmark the analysis to field data. Fortunately, these data were 
available for the commercial turboprop gearbox used in this study. 

No two gearboxes are expected to operate in exactly the same manner. Flight variables include 
operating temperature and load. Small variations in operational load can result in significant changes in 
life. Hence, the accuracy of our calculations is dependent on how close the defined mission profile is to 
actual flight operation. 

The gearboxes are condition monitored and are removed from service upon the detection of a 
perceived component failure. At the time of removal, the gearboxes are functional. The removal precludes 
secondary damage; that is, the damage is limited to the failed component. 

Field data were collected for 64 new commercial turboprop gearboxes. From these field data, the 
resultant time to removal of each gearbox is presented in the Weibull plot of figure 5(a). The failure index 
was 59 out of 64. That is, 59 out of the 64 gearboxes removed from service were considered failed. For 
these data, there was not a breakdown of the cause for removal or the percent of each component that had 
failed. The resultant L\o life from the field data was 5,627 hr and the Weibull slope e was 2. 1 89. Using the 
Lundberg-Palmgren method (above), the predicted L 10 life was 774 hr and the Weibull slope e was 1 . 125. 
The field data suggest that the L w life of the gearbox was under predicted by a factor of 7.56. 

Refurbished gearboxes 

When a gearbox is removed for cause, the failed and/or damaged components are removed. The 
gearbox is refurbished and put back into service. The 59 failed gearboxes were refurbished and placed 
back into service. The lives of these second-run gearboxes are presented in the Weibull plot of 
figure 5(b). The failure index was 34 out of 59. The times shown are based on setting the clock back to 
zero when these gearboxes were placed back into service. The resultant L w life of these second-run 
gearboxes was 1,334 hr. The Weibull slope e was 1.313. A comparison of the predicted lives and the lives 
of the new and refurbished gearboxes is shown in figure 5(c). The L i 0 life of the refurbished gearbox is 
approximately 1.7 times that originally predicted for new gearboxes but 76 percent less than that of the 
first-run (new) gearboxes. 


Reevaluation of Bearing Load-Life Exponent p 

Although errors in the assumed operating profile of the gearbox may account for the difference 
between actual and predicted life, it is suggested that using the Lundberg-Palmgren equations results in a 
life prediction that is too low for the bearings. 

With reference to equation (1), in their 1952 publication (ref. 20), Lundberg and Palmgren proposed a 
load- life exponent p = 10/3 for roller bearings, where one raceway has point contact and the other 
raceway has line contact. The 10/3 load-life exponent has been incoiporated in the ANSI/ABMA/1SO 
standards first published in 1953 (refs. 1 1 to 13). Their assumption of point and line contact may have 
been correct for many types of roller bearings then in use. However, it is no longer the case for most roller 
bearings manufactured today, and most certainly not, for cylindrical roller bearings. The computer 
analysis employed for the bearing life calculations used a value of/? = 3 for ball bearings and p = 4 for 
roller bearings. Poplawski, Peters, and Zaretsky (ref. 29) suggest the use ofp = 4 for ball bearings and 
p = 5 for roller bearings. 
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Time, hr 



Figure 5, — Field data for lives of new and refurbished 
turboprop gearboxes compared with predicted lives using 
strict-series system reliability, (a) New gearbox lives. 

(b) Refurbished gearbox lives, (c) Comparison of pre- 
dicted lives and lives of new and refurbished gearboxes. 
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From equation (11), assuming that the bearing system, which is the shortest lived component in 
the gear box and has the same Weibull slope as that of the gearbox (that is e = 2.189), 


1 1 1 
- + - 


T e s T e B T e G 
^ sys 


and substituting in the known values 


1 1 
+ - 


(5,627) 2189 L| 189 (16,680) 

and solving for the bearing system life, results in a value of 


, 2.5 


L b = 5,627 hr 


From Lundberg-Palmgren (ref. 5), the predicted bearing system life is 


f C, ' 


p 

\ eqj 


11 A hr 


then. 


f C ^ 

p 

\ eq J 


■5.21 


Calculating a revised value for the load-life exponent p for the gearbox bearings based on the actual 
bearing system life of 5,627 hr, 


( 




C r 


eq ) 


P 


~ (5. 2l) p ~ 5,627 In- 


Solving for load-life exponent p. 


p = 5.2 

where according to Poplawski, Peters, and Zaretsky (ref. 30), 

4 < p < 5 

The apparent load-life exponent p for the roller bearings is equal to 5.2. Were the roller bearing lives to be 
recalculated using the Lundberg-Palmgren theory with a load-life exponent p = 5.2, the predicted L l0 life 
of the gearbox would be equal to the actual life obtained in the field, 5,627 hr. It should be noted that if an 
exponent p = 5 were used, the predicted L i0 life of the gearbox would be 4,065 hr. This result suggests a 
strong reliance of the predicted bearing life upon the load-life exponent p. The values of the load-life 
exponent p for roller bearings equal to 10/3 from the ANSI/ABMA and ISO standards (refs. 21 to 23) and 


(14a) 

(14b) 

(14c) 

(15a) 

(15b) 

(16a) 

(16b) 

(16c) 
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4 from computer codes may provide predicted roller bearing lives that are too conservative for design 
purposes. 


Component Life 

With reference to equation (11), failure trends in the percentage of individual components failed with 
respect to the total system can be derived from the Lundberg-Palmgren model for system failure 
(ref. 5) and are expressed by Johnson (ref. 19) as follows: 

T e 

v _ ~sys 


where X is the fractional percent of components failed, L sys is the system L w life, and L„ is the component 
Co life. If the percent of each component failed as a percentage of the total number of gearboxes failed is 
known, then the L l0 life of that component can be determined from equation (17) and vice versa. 
However, a condition precedent for using equation (17) is that the individual Weibull slopes must be 
known or assumed with reasonable engineering certainty. 

Based on the Monte Carlo simulation, the component failure distribution of the gearbox for 744,450 
gearboxes is summarized in table 4. These numbers are compared with the calculated failure distribution 
predicted using the predicted component L w lives obtained using strict-series system reliability. These 
numbers are nearly identical and prove the comparability of these methods. 


TABLE 4.— PREDICTED FAILURE DISTRIBUTION BASED ON MONTE CARLO STUDY AND 
STRICT-SERIES SYSTEM RELIABILITY 


Component description 
(fig- 1) 

| Components failed, percent 

Monte Carlo study 

Strict-series system reliability 
(From eq. (17) and table 1) 

| Rolling-element bearings 

Cylindrical roller bearings 




Front pinion 

2.26 

2.45 

Rear pinion 

2.42 

2.40 

Main drive 

1.93 

1.88 

Carrier support 

0.11 

0.12 

Propeller radial 

0.64 

0.65 

Propeller thrust ball bearing 

1.45 

1.71 

Planet double-row spherical roller bearing set, 5 bearings 

89.46 

90.72 

! Gears 

Pinion 

0.11 

0.00 

Ring 

0.00 



Sun 

1.21 



Main drive 

0.01 



Planet gear set, 5 gears 

0.42 

1 

r 


On the basis of the predicted percentage of components failed from the Monte Carlo analysis and 
equation (17), the lives of the individual components are predicted based on the L 10 life of the new 
gearboxes in the field, where L sys = L i 0 = 5,627 hr. For these calculations, the Weibull slope e for all the 
components was assumed to be the same as the resultant Weibull slope of the gearboxes in the field, 
where e - 2.189 (table 3). These lives are compared with those originally predicted for the components 
from table 1. Regarding the L w lives, although the differences between those predicted from 
equation (17) for the various components are within reasonable engineering limits, there are two 
exceptions: The Co life of the five-planet, double-row spherical roller set is over seven times that 
predicted. Also, the calculated Lio life from equation (17) for the ring gear is 98.5 percent less than that 
predicted. However, the resultant ring gear life is sufficiently high enough to have only a nominal effect 
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on the gearbox life. The increase in L\q life of the five-planet, double-row spherical roller set is of 
sufficient magnitude to significantly increase the overall gearbox life from the predicted value of 774 hr 
(table 1) to 5,627 hr (table 3). 


Summary of Results 

Two computational models to determine the fatigue life and reliability of a commercial turboprop 
gearbox are compared with each other and with field data. These models are (1) the Monte Carlo 
simulation of randomly selected lives of individual bearings and gears comprising the system and (2) the 
two-parameter Weibull distribution function for bearings and gears comprising the system using strict- 
series system reliability to combine the calculated individual component lives in the gearbox. The Monte 
Carlo simulation consisted of the virtual testing of 744,450 gearboxes. These results were compared with 
each other and with two sets of field data obtained from 64 gearboxes that were first-run to removal for 
cause, were refurbished, were placed back in service, and were second-run until removal for cause. A 
series of equations was empirically developed from the Monte Carlo simulation to determine the 
statistical variation in predicted life and Weibull slope as a function of the number of gearboxes failed. 

The following results were obtained: 

1 . The resultant L 10 life from the field data was 5,627 hr and the Weibull slope e was 2. 1 89. Using 
the Lundberg-Palmgren theory, the predicted L l0 life was 11 A hr and the Weibull slope e was 

1. 125. The field data suggest that the life of the gearbox was under predicted by a factor of 7.56. 

2. The resultant L 10 life of the second-run (refurbished) gearboxes was 1,334 hr. The Weibull slope e 
was 1.313. The L\ (] life of the refurbished gearbox is approximately 1.7 times that originally 
predicted for new gearboxes but is 76 percent less than that of the first-run (new) gearboxes. 

3. From the Monte Carlo simulation, the median values for the L l0 gearbox lives equal 757 hr and 
reasonably correlate with the predicted value using strict-series system reliability of 774 hr. Flalf 
of the gearbox L w lives are less than this value and the other half more. The variation in the 
resultant life from that predicted decreases as the number of failed gearboxes in the randomly 
selected group increases. 

4. The apparent load-life exponent p for the roller bearings is equal to 5.2. Were the lives of the 
roller bearings to be recalculated using the Lundberg-Palmgren theory with a load-life exponent 
p = 5.2, the predicted L 10 life of the gearbox would be equal to the actual life obtained in the field, 
5,627 hr. If an exponent p = 5 were used, the predicted L w life of the gearbox would be 4,065 hr. 
This result suggests a strong reliance of the predicted life upon the load-life exponent p. 

5. The component failure distribution of the gearbox from the Monte Carlo simulation was nearly 
identical to that obtained by using the strict-series system reliability analysis, proving the 
comparability of these methods. 


References 

1. Zaretsky, E.V. (1992), “STLE Life Factors for Rolling Bearings,” STLE SP-34, STLE, Park Ridge, 

IL. 

2. Weibull, W. (1939), “A Statistical Theory of the Strength of Materials,” 
Ingenioersvestenskapsakademiens. Handl., 151. 

3. Weibull, W. (1939), “The Phenomenon of Rupture of Solids,” Ingeniorsvetenkapsakademiens. 
Handl., 153. 

4. Weibull, W. (1951), “A Statistical Distribution Function of Wide Applicability,” J. Appl. Mech., 
Trans. ASME, 18, 3, pp. 293-297. 

5. Lundberg, G., and Palmgren, A. (1947), “Dynamic Capacity of Rolling Bearings,” A. Polytechn., 
196. 


NASA/TM— 2007-215019 


17 



6. Coy, J.J., Townsend, D.P., and Zaretsky, E.V. (1976), “Dynamic Capacity and Surface Fatigue Life 
for Spur and Helical Gears,” J. Lubr. Technol., Trans. ASME, 98, 2, pp. 267-276. 

7. Coy, J.J., Townsend, D.P., and Zaretsky, E.V. (1975), “Analysis of Dynamic Capacity of Low- 
Contact-Ratio Spur Gears Using Lundberg-Palmgren Theory,” NASA TN D-8029, NASA, 
Cleveland, OH. 

8. Coy, J.J., and Zaretsky, E.V. (1975), “Life Analysis of Helical Gear Sets Using Lundberg-Palmgren 
Theory,” NASA TN D-8045, NASA, Cleveland, OH. 

9. Coy, J.J., Townsend, D.P., and Zaretsky, E.V. (1983), “An Update on the Life Analysis of Spur 
Gears,” in Advanced Power Transmission Technology, NASA, Cleveland, OH, NASA CP-2210, 
pp. 421-434. 

10. Savage, M., Knorr, R.J., and Coy, J.J. (1982), “Life and Reliability Models for Helicopter 
Transmissions,” NASA TM-82976 (AHS-RWP-16), NASA, Cleveland, OH. 

11. Savage, M., Paridon, C.A., and Coy, J.J. (1983), “Reliability Model for Planetary Gear Trains,” J. 
Mech. Transm. Autom. Des., 105, 3, pp. 291-297. 

12. Savage, M., Brikmanis, C., Lewicki, D.G., and Coy, J.J. (1988), “Life and Reliability Modeling of 
Bevel Gear Reductions,” J. Mech. Transm. Autom. Des., Trans. ASME, 110, 2, pp. 189-196. 

13. Savage, M., Radii, K.C., Lewicki, D.G., and Coy, J.J. (1989), “Computerized Life and Reliability 
Modeling for Turboprop Transmissions,” J. Propul. P., Trans. A1AA, 5, 5, pp. 610-614. 

14. Savage, M. (1990), “Drive System Life and Reliability,” in Rotorcraft Drivetrain Life, Safety and 
Reliability, AGARD, Neuilly-Sur-Seine, France, AGARD-R-775, pp. 35-71. 

15. Savage, M. (1991), “Life and Dynamic Capacity Modeling for Aircraft Transmissions; Final Report” 
NASA CR-4341 (AVSCOM-TR-90-C-027), NASA, Cleveland, OH. 

16. Lewicki, D.G., Black, J.D., Savage, M., and Coy, J.J. (1986), “Fatigue Life Analysis of a Turboprop 
Reduction Gearbox,” J. Mech. Trans. Autom. Des, Trans. ASME, 108, 2, pp. 255-262. 

17. Mooney, C.Z. (1997), Monte Carlo Simulation, Sage Publications, Thousand Oaks, CA, 6. 

18. Vlcek, B.L., Hendricks, R.C., and Zaretsky, E.V. (2003), “Determination of Rolling-Element Fatigue 
Life From Computer Generated Bearing Tests,” Trib. Trans., 46, 4, pp. 479-493. 

19. Johnson, L.G. (1964), The Statistical Treatment of Fatigue Experiments, Elsevier Publishing Co., 
Amsterdam. 

20. Lundberg, G., and Palmgren, A. (1951), “Dynamic Capacity of Roller Bearings,” A. Polytechn., 2, 4. 

21. (1990), “Rolling Bearings — Dynamic Load Ratings and Rating Life, First Edition,” 1S0281, 
International Organization for Standards, Geneva. 

22. (1990), “Load Ratings and Fatigue Life for Ball Bearings,” ANSI/AFBMA-9, American Bearing 
Manufacturers Association, Washington, DC. 

23. (1990), “Load Ratings and Fatigue Life for Roller Bearings,” ANSI/AFBMA-1 1, American Bearing 
Manufacturers Association, Washington, DC. 

24. Townsend, D.P., Coy, J.J., and Zaretsky, E.V. (1978), “Experimental and Analytical Load-Life 
Relation for A1S1 9310 Steel Spur Gears,” J. Mech. Des., Trans. ASME, 100, 1, pp. 54-60. 

25. Palmgren, A. (1924), “The Service Life of Ball Bearings,” NASA TT-F-13460, Translation from 
Z. Ver. Dtsch. Ing., 68, 14, pp. 339-341. 

26. Langer, B.F. (1937), “Fatigue Failure From Stress Cycles of Varying Amplitude,” J. App. Mech., 
Trans. ASME, 4, 4, pp. A160-A162. 

27. Miner, M.A. (1945), “Cumulative Damage in Fatigue,” J. Appl. Mech., Trans. ASME, 12, 3, 
pp. A159-A164. 

28. Vlcek, B. L., Hendricks, R.C., and Zaretsky, E. V. (2004), “Monte Carlo Simulation of Sudden Death 
Bearing Testing,” Trib. Trans., 47,2, pp. 188-199. 

29. Poplawski, J.V., Peters, S.M., and Zaretsky, E.V. (2001), “Effect of Roller Profile on Cylindrical 
Roller Bearing Life Prediction — Part 1: Comparison of Bearing Life Theories,” Trib. Trans., 44, 3, 
pp. 339-350. 


NASA/TM— 2007-215019 


18 



REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the 
data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this 
burden, to Department of Defense, Washington Headquarters Services, Directorate for Information Operations and Reports (0704-0188), 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302. 
Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to any penalty for failing to comply with a collection of information if it does not display a currently valid OMB 
control number. 

PLEASE DO NOT RETURN YOUR FORM TO THE ABOVE ADDRESS. 

1. REPORT DATE (DD-MM-YYYY) 2. REPORT TYPE 

01-11-2007 Technical Memorandum 

3. DATES COVERED (From - To) 

4. TITLE AND SUBTITLE 

Determination of Turboprop Reduction Gearbox System Fatigue Life and Reliability 

5a. CONTRACT NUMBER 

5b. GRANT NUMBER 

5c. PROGRAM ELEMENT NUMBER 

6. AUTHOR(S) 

Zaretsky, Erwin, V.; Lewicki, David, G.; Savage, Michael; Vlcek, Brian, L. 

5d. PROJECT NUMBER 

5e. TASK NUMBER 

5f. WORK UNIT NUMBER 

WBS 561581.02.07.03 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
John H. Glenn Research Center at Lewis Field 
Cleveland, Ohio 44135-3191 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

E- 15467-1 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 

10. SPONSORING/MONITORS 
ACRONYM(S) 

NASA 

11. SPONSORING/MONITORING 
REPORT NUMBER 

NAS A/TM-2007-2 15019 

12. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified-Unlimited 

Subject Categories: 01, 05, 37, and 39 

Available electronically at http://gltrs.grc.nasa.gov 

This publication is available from the NASA Center for AeroSpace Information, 301-621-0390 

13. SUPPLEMENTARY NOTES 

14. ABSTRACT 

Two computational models to determine the fatigue life and reliability of a commercial turboprop gearbox are compared with each other and 
with field data. These models are (1) Monte Carlo simulation of randomly selected lives of individual bearings and gears comprising the 
system and (2) two-parameter Weibull distribution function for bearings and gears comprising the system using strict-series system 
reliability to combine the calculated individual component lives in the gearbox. The Monte Carlo simulation included the virtual testing of 
744,450 gearboxes. Two sets of field data were obtained from 64 gearboxes that were first-run to removal for cause, were refurbished and 
placed back in service, and then were second-run until removal for cause. A series of equations were empirically developed from the Monte 
Carlo simulation to determine the statistical variation in predicted life and Weibull slope as a function of the number of gearboxes failed. 

The resultant L w life from the field data was 5,627 hr. From strict-series system reliability, the predicted Z, 10 life was 774 hr. From the 
Monte Carlo simulation, the median value for the L 10 gearbox lives equaled 757 hr. Half of the gearbox L 10 lives will be less than this value 
and the other half more. The resultant L m life of the second-run (refurbished) gearboxes was 1,334 hr. The apparent load-life exponent p for 
the roller bearings is 5.2. Were the bearing lives to be recalculated with a load-life exponent p equal to 5.2, the predicted L 10 life of the 

gearbox would be equal to the actual life obtained in the field. The component failure distribution of the gearbox from the Monte Carlo 
simulation was nearly identical to that using the strict-series system reliability analysis, proving the compatibility of these methods. 

15. SUBJECT TERMS 


Reliability; Life prediction-bearings; Gears; Gearbox; Probabilistic life prediction; Rolling-element bearings; Turboprop gearbox 


16. SECURITY CLASSIFICATION OF: 

17. LIMITATION OF 
ABSTRACT 

18. NUMBER 
OF 

19a. NAME OF RESPONSIBLE PERSON 

STI Help Desk (email:help@sti.nasa.gov) 

a. REPORT 

u 

b. ABSTRACT 

u 

c. THIS 
PAGE 

u 

uu 

PAGES 

24 

19b. TELEPHONE NUMBER (include area code) 
301-621-0390 


Standard Form 298 (Rev. 8-98) 
Prescribed by ANSI Std. Z39-18 
































